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We report resonant inelastic x-ray measurements of insulating La2Cu04 and Sr2Cu02Cl2 taken 
with the incident energy tuned near the Cu K absorption edge. We show that the spectra are well 
described in a shakeup picture in 3rd order perturbation theory which exhibits both incoming and 
outgoing resonances, and demonstrate how to extract a spectral function from the raw data. We 
conclude by showing q-dependent measurements of the charge transfer gap. 
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Inelastic x-ray scattering (IXS) has shown promise as 
a practical probe of electronic excitations in condensed 
matter because of its broad kinematic range and direct 
coupling to the electron charge. However, since x-rays are 
strongly absorbed in high density materials, successful 
applications of the technique have been limited to low-Z 
systems 

Several recent studies @|j]||f§, have shown that by 
that by tuning the incident energy near an x-ray absorp- 
tion edge a Raman effect could be measured, despite the 
high absorption, because of the resonant enhancement. 
These studies have emphasized the role of coulomb inter- 
actions in the scattering process. Since it involves cou- 
pling between highly excited virtual states and strongly 
correlated valence states, it is important to characterize 
the resonance process well for the technique to be useful. 

With emphasis on systematics, we have measured res- 
onant inelastic x-ray scattering at moderate resolution 
(Ai5=0.9 eV) near the CuK absorption edge in the high- 
T c parent insulator La2Cu04 (LCO) as a function of in- 
cident photon energy. Based on the changes of inelas- 
tic intensity and peak position with incident energy we 
show that the scattering is well described in a "shakeup" 
picture in 3rd order perturbation threoy |Tc|| . We also 
present higher resolution measurements (Ai?=0.45 eV) 
on another insulator, S^CuC^C^ (SCOC), as a function 
of momentum transfer, q, which show some new features, 
such as the 2 eV optical charge transfer gap. 

Experiments were carried out at the X21 wiggler line at 
the National Synchrotron Light Source and the 3ID (SRI- 
CAT) undulator line at the Advanced Photon Source. 
At X21 the energy resolution was 0.45 eV and typical 
count rates were 0.4 Hz. At 3ID with 0.9 eV resolution 
9 Hz was typical. The scattered light was collected with 
a spherical, diced, Ge(733) analyzer and imaged onto a 
detector. Energy analysis was done by rotating the an- 
alyzer and translating the detector in coincidence. The 
momentum transfer was varied by rotating the entire ap- 
paratus around the scattering center (exact experimental 



geometries are indicated in the figure captions). 
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FIG. 1. Resonance Profile for La2Cu04 taken at 3ID with 
e.i\\a and q=1.27A _1 parallel to c. The upper frame shows the 
raw spectra plotted against transferred energy (curves are off- 
set for clarity). In the lower frame the open and filled circles 
are the inelastic peak height and position, respectively, plot- 
ted against incident energy. The black line is the fluorescence 
yield which peaks at the Culs — > 4p energy. 



The LCO and SCOC crystals were grown by techniques 
described previously They were characterized 

with a spectroscopic ellipsometer to assure surface qual- 



1 



ity. 

In the upper panel of Fig. 1 we show inelastic x-ray 
spectra from La2Cu04 at fixed q for 16 different incident 
x-ray energies across the Cu K absorption edge. Since the 
overall absorption change across the edge is only 14% (x- 
ray absorption being dominated by the La atoms) cor- 
rections for sample self-absorption did not significantly 
alter the line shapes. So shown here are the raw spectra. 
The strong peak at zero energy loss is elastic scattering. 
The energy gain side (negative energy loss) shows a back- 
ground of 12 counts per minute. The principal feature of 
each curve is a single peak, seen previously in Nd2CuC>4 
H, whose position, intensity, and line shape vary greatly 
as a function of incident energy. 

We summarize the resonant behavior in the lower 
frame of Fig. 1 where the inelastic peak height (open 
circles) and position (filled circles) are plotted against 
incident energy. The thick line is the CuK fluorescence 
yield which shows the location of the edge (a localized 
Culs — > 4p transition). The peak height shows two max- 
ima, the stronger of which correlates with the peak of the 
edge and the weaker with the pre-edge shoulder. In both 
cases the maximum is offset from the absorption peak by 
about 2.5 eV. 

The peak position shifts nonlinearly with incident en- 
ergy; below and above a resonance it is roughly linear 
while near a resonance it plateaus. This behavior differs 
fundamentally from a classic Raman effect, in which one 
expects a linear peak shift with unit slope with a gradual 
rise in intensity below a resonance, and saturation above 

Independent of any model, this type of low energy loss 
scattering leads to excited states of the valence electron 
system in the absence of core excitations. However, be- 
cause one is near the Is absorption threshold of Cu the 
scattering proceeds through a set of almost real, highly 
excited (9 keV) intermediate states, which have a Is core 
hole and an extra electron excited in a localized Ap state 
(IsAp). When the highly excited intermediate state dis- 
appears it leaves behind low-lying valence excited states 
- in principle conserving energy and crystal momentum. 

One must describe these many-body intermediate 
states, i.e. their matrix elements as well as their off-shell 
weight, in order to characterize the coherent, second or- 
der process. Different groups have resorted to different 
approximation schemes. Starting with N valence elec- 
trons in a small cluster, Tanaka and Kotani Jl4|] describe 
the intermediate states as a set of N+l interacting elec- 
trons in the presence of a rigid impurity - the Anderson 
Impurity Model. This treatment assumes that the core 
state is suddenly created, and that it can be treated as 
a fixed, external potential. It emphasizes the multiplct 
coupling among the N+l electrons, and is done numer- 
ically in exact diagonalization with a large number of 
basis functions. 

Taking a more analytic approach, Platzman and Isaacs 



[ |i0[ treat the many-body problem by describing Coulomb 
interactions among electrons in a series of perturbation 
diagrams p5| . This approach assumes that interactions 
can be taken to be weak for a suitable choice of basis 
functions or that one can sum enough terms in perturba- 
tion theory to include the important physics. They argue 
that near a sharp, dipole-allowed transition the dominant 
term occurs in 3rd order. Writing it out explicitly one 
arrives at 
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In this expression the sum is on all possible states of the 
Is hole and 4p electron, Ef s4p is their energy, and jk 
is their inverse lifetime. The numerator contains matrix 
elements for absorption, M a i, s — (e/mc)(ls4p|p • A|k;}, 
emission, M em = (e/mc)(k s |p ■ A|ls'4p'}, and coulomb 
interaction between the core and valence states, M coll / = 
/ dxdx' {Is 4p'; /|p(x)/S(x')/|x - x'||ls4p;i). kj and k s 
are the incident and scattered photon momenta, (i.e. 
q=ki-k s ) and uj = u>i — uj s is the energy loss. 

Physically this expression represents the following. 
The incident photon, with energy tuned to the CuK ab- 
sorption edge, creates a virtual IsAp pair on a Cu site. 
This pair is bound as an exciton by the coulomb interac- 
tion (not included in Eq. (1)) and so is non-dispersive. 
It takes up the momentum of the incident photon, 4.55 
A -1 , and scatters off the valence electron system. When 
the exciton recombines, the emitted photon reflects the 
energy and momentum imparted to the system. This 
is commonly called a "shakeup" process, which to first 
order in the coulomb interaction is given by Eq. (1). 

To get a transition rate we square the quantity (1) and 
perform an incoherent sum on final states. We postulate 
that the intermediate states are approximately degener- 
ate with energy Ef s4:p = Ek (since they are spatially 
localized), which allows factoring of the energy denomi- 
nators from the sum. We arrive at 
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and 6i and e s are the polarizations of the incident and 
scattered photons. The two lorentzians in (2) are incom- 
ing and outgoing resonances in the photon frequency, so 
we see that this treatment is completely analogous to 
third-order optical Raman scattering from phonons in 
semiconductors, in which the scattering is described by 
a single operation of the electron-phonon interaction on 
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a virtual valence electron- hole pair |lq|. The two reso- 
nances formally come about the same way. 

The central result of this paper is Eq. (2). It says that, 
within our approximation, all the very different spectra 
in Figure 1 derive from the same fundamental quantity, 
Sk, which depends on the difference u> = LUi — uj s rather 
than on Ui and lo s independently. A way to test this 
result would be to take the curves from Fig. 1, divide 
each by its respective denominator from Eq. (2), and see 
if they all collapse to the same function. 
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FIG. 2. The result of dividing energy denominators 
in (2) from the experimental spectra in Fig. 1. For 
E K = (8995.14 ± 1.61) eV and lK = (2.38 ± 0.542) eV all 
spectra collapse onto a single curve. The thick line is a Gaus- 
sian fit. 

In Eq. (2) we assumed that we are near a single res- 
onance, so we take the nine highest curves from Fig. 1 
(around the second peak in the resonance profile) and 
subtract their background and elastic scattering. We use 
Ek and jk as flexible parameters and divide each spec- 
trum by its respective denominator. Using a nonlinear 
fitting algorithm, we adjust the values of Ek and "fx 
to minimize the total variation (the x 2 summed over all 
points and all spectra) among the curves, irrespective of 
the resulting shape. For the values Ek = (8995.14±1.61) 
eV and jk — (2.38 ± 0.542) eV we find that the spectra 
collapse onto a single curve, shown in Figure 2. 

The result for Sk (q, w; e s ) is a single peak at 6.1 eV 
energy loss and width of 3.9 eV. Referring to the cluster 
calculations of Simon |l7]] we suggest that this feature 
is a transition from the b\ g ground state to an a\ g exci- 
tonic excited state composed of symmetric combinations 
of a central Cu3c? a ,2_ J/ 2 orbital and the surrounding 02p a 
orbitals. 

To illustrate what qualitative aspects of the data are 
captured by our fit, i.e. by the resonant denominators 
in (2), we take a single function for Sk, i-e. a fit to the 
collapsed data in Fig. 2 (thick line), combine it with 
the denominators in Eq. (2), and produce the model 



resonance profile shown in Fig. 3 (identical formatting 
to Fig. 1). The salient features are reproduced, including 
the peak shift with incident energy and the 2.5 eV offset. 
All this behavior comes from the energy denominators in 
(2) and is independent of the nature of the core state or 
the particular valence excitation. 
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FIG. 3. A model resonance profile generated from the dark 
line in Fig. 2 and the two resonant denominators in equation 
(2). Formatting is identical to Fig. 1. The solid line is a sim- 
ulated CuK fluorescence yield showing the absorption peak. 
The continuum above the peak is omitted since it does not 
contribute to the resonance. 

In this simple shakeup description Sif(q, w; <ii, e s ) has 
an explicit relationship with S{(\ 1 lu), the dynamical 
structure factor measured in nonresonant inelastic x-ray 
scattering jl8j . This can be seen by writing out the ma- 
trix element M cou i in momentum space, which (neglect- 
ing exchange between core and valence states) has the 
form 

M cou i = V 4 ^ 2 -Ff 84? ,(q + G; ei)(f\p v ^ +G \i). (4) 

Here Fi s4p (k) is the static x-ray form factor of the ls4p 
exciton. It is dependent implicitly on the incident polar- 
ization ej since in the dipole approximation M a b s deter- 
mines the spacial orientation of the Ap. p v ^ is the valence 
part of the many body density operator p q and the sum 
in (4) is on all reciprocal lattice vectors, G. 

The quantity (f\p v ,q\i), when squared, multiplied by 
5(u> — Ef + Ei), and summed on final states, is iden- 
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tically the valence part of the dynamic structure factor 
S(q, uj) . Doing the sum on G before squaring we find that 
Si<-(q, u>; ej, e a ) is a superposition of many 5(q + G,w) 
functions, weighted by the form factor of the core states. 
Therefore, Sk can be thought of as a response function 
similar to S "projected" onto the form factor of the core 
state. 
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FIG. 4. Resonant spectra from Sr2Cu02Cl2 as a function 
of q, taken at X21 with both q and E approximately parallel 
to the (210) crystallographic direction. The feature at 2 eV is 
the charge transfer gap, which changes shape and shifts with 
increasing q. 



Finally we present some higher-resolution (AE=0.45 
eV), q-dependent spectra from the insulator Si^CuC^C^ 
(Fig. 4). The salient spectral features are three peaks: 
the same feature we saw in LCO (appearing here at about 
5 eV), a feature at 2 eV which appears to shift and lose 
definition with q, and a peak at 10.5 eV which is ab- 
sent at low-q but which gains intensity as q is increased. 
The 2 eV feature strongly resembles dipole-active charge 
transfer gap seen with optical reflectivity Jlj| . 

To summarize, we measured RIXS spectra at the CuK 
edge in LCO and SCOC. From the resonance profile we 
deduce that the scattering can be described as a shakeup 
process in 3rd order which, analogous to optical Raman 
scattering from phonons, exhibits incoming and outgo- 
ing resonances. The scattered intensity has the form of 
a doubly resonant form factor multiplied by a response 
function, Sif(q, w; ti, e s ), which can be thought of as the 
dynamical structure factor projected onto the form factor 
of the intermediate core state. 

The advantages of this description are (i) that it makes 
no assumption that the core state is rigid and so mo- 
mentum conservation enters naturally, and (ii) that it 
allows one, given certain knowledge of the intermediate 
state (i.e. Ff s4p (k)), to relate the scattering to a response 
function in terms of the valence electrons only. This de- 
scription is useful when the core resonance is sharp and 



well isolated, and when one is mostly interested in the 
valence electron spectrum and is willing to sacrifice a de- 
tailed multiplet description of the core states. 
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